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Abstract 

The role of spin-orbit interaction, arises from the Dzyaloshinski-Moriya anisotropic antisymmet- 
ric interaction, on the entanglement transfer via an antiferromagnetic XXZ Heisenberg chain is 
investigated. From symmetrical point of view, the XXZ Hamiltonian with Dzyaloshinski-Moriya 
interaction can be replaced by a modified XXZ Hamiltonian which is defined by a new exchange 
coupling constant and rotated Pauli operators. The modified coupling constant and the angle of 
rotations are depend on the strength of Dzyaloshinski-Moriya interaction. In this paper we study 
the dynamical behavior of the entanglement propagation through a system which is consist of a 
pair of maximally entangled spins coupled to one end of the chain. The calculations are performed 
for the ground state and the thermal state of the chain, separately. In both cases the presence of 
this anisotropic interaction make our channel more efficient, such that the speed of transmission 
and the amount of the entanglement are improved as this interaction is switched on. We show that 
for large values of the strength of this interaction a large family of XXZ chains becomes efficient 
quantum channels, for whole values of an isotropy parameter in the region — 2 < A < 2. 
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I. INTRODUCTION 



Recently transmitting a quantum state is a most important task in quantum information 
and computation processing Today this purpose could be achieved in two manner: i) 
using standard teleportation protocols and, ii) information transfer via quantum networks. 
Bayat et.al [2] have shown that the fidelity of transmission is the same for both cases. 
Additionally, the former case employs the flying qbits as quantum channel and hence the 
fidelity of transmission reduces due to incompleteness of stationary-to-flying qbits conversion 
process [3|. Thus the later case is superior to the former, particularly in solid state devices. 
Among the numerous quantum systems suitable for quantum networks implementation, the 
spin chains offer a great advantages. One of the most interesting art of the spin chains is 
their ability to use as quantum wires in the information transfer protocols over the short 



distances 



1 If ] . Tunable spin interaction in these systems plays the key role to motivate 



one for using this permanently potential in the quantum information transfer processing. 
The major works on spin chain are in the Ferromagnetic(FM) phase [7| and the effects of 
temperature [8j and decoherence j{| have been investigated for FM channels. There are 
lesser works on Antiferromagnetic(AFM) phase 12j, ll3(, but as mentioned in Refs. 
AFM spin chain have higher ability in transfer of information and thus is a better and 
faster alternative. So we prefer to study AFM spin chain. Fortunately, Antiferromagnetic 
spin chains with short length(up to 10 spins) have been built experimentally 14|. However, 
much attention has been paid to the entanglement in spin systems with only spin-spin 
interaction and spin-orbit interaction has been leaved. The spin-orbit interaction produces 
and anisotropic part of exchange interaction between localized conduction band electrons 
in crystals with lack of inversion symmetry, including all low dimensional structures and 
also bulk semiconductors with zinc-blende and wurtzite type of crystal lattice 15j. The 
main part of the anisotropic interaction has the form of Dzyaloshinski-Moriya interaction 
(DMI) [l6-18] which have explained the weak ferromagnetism of antiferromagnetic crystals 
(a — Fe^Oj,^ MnCO% and CrF%). As was shown 19l-l2lj in the two-qubit Heisenberg systems 
the (DMI) plays an important role. The (DMI) express by 

This interaction arises from extending the Anderson's theory of superexchange interaction 



by including the spin-orbit coupling effects 



17(. In this paper we interested to consider the 



model with DMI and also investigated the effects of this interaction on the information 
transfer processing. 

Following the approach proposed by S. Bosej^], we place a spin encoding the state at one 
end of the chain (which is now equipped with DMI) and wait for specific amount of time to 
let this state propagate to the other end. Therefore, entanglement could be transferred from 
one end of the chain to the other. We will show that the XXZ chain with DMI can be reduced 
to the modified XXZ chain with new coupling constant and rotated Pauli operators. The 
modified coupling constant and the angle of rotations are depend on the strength of DMI. 
Consequently, the entanglement transfer protocol through the XXZ+DMI chains becomes 
the same as a protocol via a general form of XXZ chain with new definition of coupling 
constant and Pauli matrices, but it should be noted that at this case the initial state may 
be changed. Indeed, both modified coupling constant and the angle of rotations, which is 
referred by "phase factor" through the text, are the efficient parameters for entanglement 
transfer. For the sake of clarity, the effects of phase factor and coupling constant have been 
separated. Our calculation shows that this phase factor does not affect on the entanglement 
transfer when A = for the XXZ chain. Whereas, this phase factor has desirable effects on 
the parameters of entanglement transfer for the case of A < and undesirable effects when 
A > 0. However, there is a contest between the phase factor and coupling constant. For 
large amounts of the strength of DMI(D), the effects of coupling constant are more dominant 
than the effects of phase factor and has desirable effects. Also, in the whole range of A, the 
effects of the modified coupling constants and the phase factor simultaneously investigated 
on our main goal, i.e, information or entanglement transmission. Furthermore, the other 
advantage of DMI is that the speed of information transmission increases as increasing the 
strength of this interaction. Also, at nonzero temperature this interaction (DMI) improves 
both E max and t opt as increasing the strength of DMI (D) for both positive and negative 
amounts of A. 

The paper is organized as follows. In Sec. II we introduce the XXZ Hamiltonian with DMI 
and discuss the way to change the form of Hamiltonian to the general XXZ Hamiltonian. 
The new set of eigenstates have been introduced in this Sec. In Sec. Ill we introduce our 
state transmission. In Sec. IV we show analytically the dependence of concurrence on D 
and pase factor. In Sec. V we show our numerical calculation of entanglement transfer 
and corresponding times in both zero and nonzero temperature and calculated the speed of 
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information transfer while, we summarize our results in Sec. VI. 



II. HAMILTON! AN AND MODEL 

The Hamiltonian of the open XXZ chain in presence of spin-orbit interaction is defined 

as 

N ch -1 N ch -1 

H ch = i J i^i+i + °MhJ + Aofo^J + x *»M> C 1 ) 

i=l i=l 

Where N c h is the number of spin in ID chain, Oi={of ,of ,of) is the vector of pauli matrices, 
J is the exchange coupling, the parameter A is the anisotropy exchange coupling in the z 
direction and D is the Dzyaloshinski-Moriya vector. Different phases of chain is depending 
on different range of J and A. For the case of J < the chain is called Ferromagnetic(FM) 
Heisenberg chain and the case J > is Antiferromagnetic(AFM) Heisenberg chain. The 
AFM Heisenberg chain includes FM phase(y < 0), AFM phase(0 < j- < 1) and Neel 



phase(4 > 1) 22]. If we take D = Dz, then the above Hamiltonian can be written as 



Nch-1 N ch -1 

H ch = £ {JWM + i + + MtoUi) + dJ2 K<i - ^f+J- ( 2 ) 

i=l 1=1 

N ch 

This Hamiltonian is invariant under z-axis rotation, i.e, [H,S Z ] = 0, where S z = h^of. 
By this property the Hamiltonian can be express in the form of usual XXZ model without 
explicit DMI. For this purpose, the pauli operators iu x-y directions are uranipulated by the 
unitary transformation which is depend on the spin sites [2c 

N ch 

U Nch = exp »=a , (3) 

where <fi = tang -1 (D/ J). Hence these spin coordinate transformations reads 

of = of cos 4>i + of sin </>;, 
of = —of sin 4>i + of cos i; 

5? = (4) 
where <fii = (i — 1)0. So, the modified Hamiltonian is express as 

N ch -l N ch -1 

H ch = U Nch H ch Ul h = J £ K^+i + nn + x) + A ^ a*a* +1 , (5) 

i=i i=i 

J = S gn(J)VJ 2 + D 2 . 



The eigenstates of the new Hamiltonian (|V> n )) are related to the earlier one (\ip n )) via 
\4>n) = U]y c h\4>n)- As the model of the system is specified, we can investigate the information 
transmission processing in this system. 



III. ENTANGLEMENT TRANSMISSION 



The quantum information transmission of one part of a two-spin maximally entangled 
state (O'O) via XXZ+DMI spin chain is investigated while the spin chain is in its ground state 
(\if) g8 )ch)- At t=0 we interact the spin with the first spin of the chain. We suppose that 
the chain is prepare in a unique grand state, initially. The preparation could be performed 
by applying a small magnetic field, if it is required. Furthermore, the interaction between 
spin and first spin of the chain has the same form of the rest of interaction, 



Hj = J(o%of + a»a» + A<r V*) + D(a x a\ - a y af) 



(6) 



Indeed, the system is consist of O'O and N c h spins and hence the total length of the system 
is N = N c h + 2. The initial state of the system is 



|^(0)) = |V> )o'0 ® \lf>ga)ch, 



(7) 



where 



)o'o 



|01> — |10> 
V2 ' 



(8) 



This \ip(0)) is used as a channel which transfer the entanglement. Therefore, the total 
Hamiltonian being 

H = I (y ®(H ch + H I ), (9) 

By this Hamiltonian, the initial state evolves to the state \ip(t)) = e~ lHt \i()(0)) and the two 
sites reduced density matrix of can be computed as p mn (t) = t^rnn{\4'(t)){'ip(t)\}, where tr^ 
is the partial trace over the system except sites m and n. The two sites reduced density 
matrix in computational basis ( 1 00) , 1 01) , |10), |11)) has the general form asfl^ 

/ a(t) \ 



Pmn (t) 



x(t) z(t) 
z*(t) y(t) 
\ b(t) j 



(10) 



Although, the density matrix has been written in the Schrodinger picture, but p mn (t) in 
terms of spin-spin correlation function could be expressed in the Heisenberg picture as 

follows y 



o(t) = i + i«(t) + <(t)) + «(tK(t)), 

y(t) = 1 - - <(t)> - «(*)<(*)>, (11) 

b(t) = 1 - i«(t) + <(t)> + «(*)<(*)>, 

*(*) = «W<W> + (^(i)<W} + »(KW<W}-R(*KW))- 

these correlations are computed in terms of the initial state (EqJTj) and cx^(t) = e~ %m a c ^ l e~ %m 
where a = {x, y, z}. Since the concurrence is directly defined in terms of the density matrix 
and so any minimization procedure is not necessary, the concurrence is used as a measure 



of entanglement for arbitrary mixed state of two qubits 



25|, 



C = max{0,2X max - trVM}, (12) 
R = pa y <g) a y p*a y <g> a y } (13) 



where A max is the largest eigenvalues of the matrix y R. For our density matrix the concur- 
rence results to be 

C = 2mas{0,C (1) ,C (2) }, (14) 
where = —yjxy and = \z\ — y r ab. Because, is always negative here the 



concurrence is 



22 1 



C = 2max{0, \z\ - Vd)}. (15) 

For our main goal the first site refers to 0' and other site refers to the spin located at the 
end of the chain, say (j). So in terms of density matrix the subscript m is change to 0' and 
n is change to j. The singlet fraction of the state, po'j, as an indicator of the average fidelity 
of state transferring could be obtain as 

F=(il>-\po>3\il>-) = l(x + y-2z). (16) 
With the aid of H c h = Ujsr ch H c hUlf we have 

«(*)> = W) 

= (0'0|(^ s | c ,^ ch+1 e-^<e^^ ch+1 |^) c ,|0'0), (17) 



with use of the \ip n ) = UNchl^n), the expectation value can be express as 

«(*)> = (0'0|(4 s | ch ^ c;i ^ ch+1 <(t)^ ch+1 ^ c J^ s ) c ,|0'0), (18) 

N ch 

-i £ (2m-l)iK, 

where UN ch UN ch +i = exp m=1 is phase factor which modify the states on the right 

hand of above equation. So, we can conclude that this model is similar to usual XXZ model 
with the new strength coupling in XY direction (J) and the new set of states which are 
multiplied by the phase factor. In the following, the effects of these parameters (phase 
factor and J) will investigated on entanglement transfer processing. 



IV. ANALYTICAL CALCULATION 



To more clarifying, the concurrence between the 0' site and the end of the chain with the 
length oi(N c h = 2) have been calculated analytically in the appendix. In these calculations 
the explicit form of p mn (t) = tr mn {\^{t)){^{t)\) was used and these results are quality 
compatible with the numerical results for higher N. From the relation (jXSJ) the concurrence 
for the case of (A = 0) is 

C(t) = ^(4|cos(£t)-l|-|e 2 ^(l + cos(£t)) 

o 

+ iV2 sin (ft) | x |l + cos (£t) + iV2e 2i *sm (19) 

where f = 2\/2J and we have the maximum entanglement (C=l) for ft = n. As we can see, 
amounts of entanglement at the first peak (E max ) and corresponding time (t = t opt ) in this 
case are independent on the <fi. 

Furthermore, for the case of A ^ the concurrence is obtained in Eq.f lAlOl) . In this case 
E max and t op t depend on as indicated in Fig. 1. In these figure, for clarifying the role of <p, 
individually, we fixed the size of J. The results show that, in the case of A < 0, increasing 
(f) improves both E max and t opt and for the case A > 0, the increase of has undesirable 
effects on both E max and t opt . The calculation becomes more involved when N exceeds 4, 
this prevents one from writing an analytical expression for the concurrence. So, we solve 
this problem numerically. 
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V. NUMERICAL CALCULATION 



A. Entanglement at T=0 

Whereas, in the XXZ spin chain the AFM Heisenberg chains ( J > 0) is a better candidate 
for information transfer than FM chains (J < 0) 11], we confined our calculations to the 
case J > and take J = 1 to simplify the calculation. The entanglement is calculated 
between 0' and spin located at the end of the chain which has the length of N=8. In Fig. 
2(a), E max , as measured by concurrence, have been plotted in the domain of (—2 < A < 2) 
for different values of D and Fig. 2(b) illustrates the behavior of E max as a function of A for 
different values of 0, where J is fixed. For the special case D=0 (i.e, J = J and = 0) the 
results are the same as in Ref. |2|, qualitatively. In this case E max vanishes at the quantum 
phase transition (QPT) point (4 = —1) 22]. In the presence of D, J ^ J and hence the 



QPT point shifts to the ^ = — 1, e.g, for the case of D = 1 the QPT occurs at A = — v2 



281 ] . Despite to the results of [2], at this modified transition point E max does not vanish, 
this is due to the presence of the phase factor as indicated in Fig. 2(b). Furthermore, 
the amount of J is so large for large values of D and so the QPT point disappear in the 
frame of — 2 < A < 2. In this range of A and D, the amount of E max is unsensible to the 
values of and hence J plays the main role to quantifying E max . Since, j approaches zero 
as D becomes large, the system treat as the case A = 0, i.e, E max — > E max (A = 0). Fig. 
3(a) reveals that the behavior of t opt is compatible with the result of Fig. 2(a). This figure 
shows that t opt decreases as increasing D, hence the presence of DMI enhances the speed of 
information transmission. Also, Figs. 2(b) and 3(b) show that the effect of on E max and 
t opt is undesirable for the case of A > 0. In contrast, for the case A < 0, the maximum 
entanglement at the first peak and speed of transmission enhances with 0. In summary, all 
of the chains with —2 < A < 2 can be used as protocol for information transfer processing 
with the same cost. 

Fig. 4 depicts the singlet fraction in terms of A for different values of D. The result of this 
picture confirms the above mentioned effects of DMI on entanglement transfer properties. 

In order to better illustrate the effects of D on the entanglement transfer properties, E max 
is plotted in terms of D and in Figs. 5(a)-5(d) for different values of A respectively. The 
dips appearing in Figs. 5(a) and 5(c) are due to the effects of the phase factor (0) which 



S 



is obvious from Figs. 5(b) and 5(d). The effects of DMI on speed of state transmission are 
plotted in Figs. 6(a) and 6(b), more obviously. 



B. Thermal Entanglement 

Since, preparing the system at T = is far from access, the presence of thermal excitations 
is unavoidable. So, we consider the state of channel as a thermal state instead of ground 
state. The thermal state of channel at temperature T is given by density matrix p c h = 6 ? " oh , 
where (3 = j^f, &b is the Boltzman constant and Z = tr(e~ l3Hch ) is the partition function. 
So the initial state of system is 

p(0) = |^-)(V-|®— (20) 

Employing the system parameters as before, we repeat the calculations for this new initial 
state and the results have been shown in Figs 7 and 8 for A = 1 and A = — 1, respectively. 
Two considerations are in order at this stage. Firstly, at hight temperatures, thermal fluctu- 
ations suppress the quantum correlations and hence the entanglement vanishes and also, t opt 
becomes so large. Secondly, the presence of DMI amplify the quantum correlations, so by 
increasing D we can obtain nonzero entanglement at larger temperatures. For instance Fig 
8(a) shows that for the case of A = — 1 the entanglement can be exist at higher temperatures 
in the presence of DMI while it is zero for all temperature in the absence of DMI. Also, the 
speed of transfer improves as D increases. 



C. Information speed 

As mentioned before, the DMI imposes desirable effects on the speed of transmission. In 
order to better clarify, we compare t opt with v _1 (which v is the spin-wave velocity) which is 



obtained with the aid of field theoretic techniques [27J. Note that, the later get the qualitative 
behavior of correlations in the thermodynamic limit while the former is calculated for very 
short chain (N = 8). Also, in the field theoretic techniques the dynamical correlations are 
computed for the ground state of the system while in our problem it is not the case. Despite 
difference, we still can use some well-known results of the field theoretic techniques. For 
instance, in our modified XXZ model(XXZ model equipped with DMI) and for the range 
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— 1 < A < 1, the spin-spin correlation function in the asymptotic thermodynamical limit 



have the following form 



27] 



WK+n >~ (-l)V-^]- (1/a ^, (21) 
% = % = ^ = 1 - -, 

7T 

where 

A = cos( 7 ), (22) 
and the propagation velocity of excitation in the chain can be written as 

^ ?rJsin(7) . (23) 

7 

This velocity is proportional to the strength of DMI (D), which is included in J, and also it 
depends on the size of A. According to the Ref. 27j], v refers to the velocity of propagation 
of the correlations which is related to the entanglement transmission speed. Fig. 9 shows 
V 1 , which is obtained from Eq. (123]) . and t opt , which is calculated in previous section, in 
terms of D for two different values of A. As we can see both quantities have similar behavior 
qualitatively, such that they descent with increasing of D. As this figure illustrates, despite 
of the field theory calculation for spin-wave propagation time (u _1 ), the curves of t pt cross 
each other, indeed at the cross point the effect of the phase factor (0) becomes considerable. 
As a consequence, for large values of D the speed of information transmission raises up and 
ultimately reaches the asymptotic value. Faster dynamics in large values of D stems to the 
entanglement enhancement of the channel in this region. 



VI. CONCLUSIONS 



In this paper we examined a XXZ Heisenberg chain equipped with the Dzyaloshinski- 
Moriya interaction (DMI) as a quantum channel for investigation of the entanglement trans- 
fer properties. We had shown that the presence of DMI enhances the amount of coupling 
constant, J, and also imposes a phase factor on the state of the channel. In order to clarifying 
the role of DMI, we trace the effects of the new coupling constant (J) and phase factor (0) 
for a wide range of the anisotropy parameter, A, separately. Indeed, increasing J with D 
leads to increase the strength of spin-spin correlations and ultimately improves the amount 
of the entanglement. For the case A < the effects of <fi on the entanglement properties of 
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transmission is more desirable. In contrast, for the case A > increasing 0, individually, 
destroys the entanglement of the channel. The effects of J dominate for large values of D 
and hence the channel becomes more efficient for all values of anisotropy parameter in the 
region —2 < A < 2. We calculated the entanglement properties for the ground state and the 
thermal state of the chain, separately. Our results show that the amount of entanglement 
and the speed of transmission increase as D increases. Also, we show that the entanglement 
can be exist at higher temperatures as D increases and hence the transmission channel could 
be work more efficiently at higher temperatures. 



Appendix A: Analytical Calculation for N=4 

In this appendix, we give analytical calculation for the XXX and XXZ chains with N = 4, 
separately. For XXX chain (A = 0), the initial state |^(0)) in the basis of the Hamiltonian 

m 



for N = 3 can be written as 

-i4> 



Mo)) 



0) ' ® 



2 + e 2 ^ 5 > + -(-v / 2 + e 2 ^ 7 >) 



+ \i)o>®(-^\fa) + \(i 



where 



y/2^*)\ik) + \(i 



y/2^*)\ih))].. 



(Al) 



fa) = 


-L(-|on) + |iio)), 


fa) = 


_L(_|001} + |100)), 


fa) = 


|()011> - a/2|101> + |110», 


fa) = 


^(|001)-v / 2|010) + |100)), 


fa) = 


i(|011) + V2|101) + |110)), 


fa) = 


i(|001) + v / 2|010) + |100)), 



(A2) 



are the relevant eigenstates of H 3 (A = 0) and the corresponding eigenvalues are E 2 = E 3 = 
0, E 5 = E % = — £ and E 7 = E s = £, here we define £ = 2y/2J. The state of system at later 
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times can be obtained as 

p —i4> p 2i<f> 1 

m)) = —[\0)v®(^=\ih) + -(v / 2 + e 2 ^)e* f U> 
+ \{-V2 + ^*)e^\ih)) + |l)o' ® (-^1^3) 
+ 1(1 + V2e^)e^ 6 ) + 1(1 - v^e^e-^lVs})]. (A3) 

The corresponding reduced density matrix, is in the form of Eq. ( jTUl) with the 

following components 



^|e 2 ^(l + cos (#)) + zv^sin (fi)| 



16 



x = — | e 2 ^(— 1 + cos (#)) + zV2sin(£t)| 



16 



y 



+ llv^cos^ + ze^sin^)! 2 , 

o 

= l\V2e 2 ^cos(£t)+ism(£t)\ 2 , (A4) 
8 

+ t^I - 1 + cos(£t) +iv / 2e 2 ^sin(^)| 2 , 



16 

b = h + cos(^) +iv / 2e 2 ^sin(^)| 2 , 
Id 1 1 

z = l(-l + cos(^)). 
Therefore the concurrence could be computed as 

C(t) = 1(4| - l + cos(ft)| - |e 2 ^(l + cos(£t)) 

o 

+ z\/2sin(£t)| x |l + cos(ft) 

+ zV2e 2 ^sin(^)|). (A5) 

Following the same procedure for XXZ chain (A 7^ 0), the initial state \i/)(0)) in the basis of 
corresponding Hamiltonian is 



|, (0) ) = _ [ |o) ,«(_ W + ^=|, 5 ) + -^ F |« ) 



+ \^®(-^s) + ^^\A) + ^==m)], (A6) 
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where a = A+ ^f^ = _ Wp^i and 



2.J ' ' 2J 

\rh) = ^=(-|011) + |110», 

M = -J=(-|ooi) + |ioo», 

^ = ^i=^(|oii)-«|ioi) + |iio)), 

V 2 + cr 

V 2 + cr 

l^r) = - 7 =L=(|011)+^|101) + |110)), 

= ^^ (I001> + /9|010) + |100)), 

are the relevant eigenstates of H 3 (A ^ 0) and the corresponding eigenvalues are E 2 = E 3 = 
0, E 5 = E 6 = —2Ja and E 7 = E 8 = 2.1(5. Therefore, the initial state, |^(0)) evolves to the 
state 

2i<j> ( n _|_ e 2i(/>^ e -iE 5 t 



(_o _i_ p 2i(f>\ p -iBht 1 
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The corresponding reduced density matrix, /?q/ 2 (£) is in the form of Eq. (jTUl) with the 
following components 



and hence 



1 e 2iJta (a + e 2i ^) 1 2i(f) e - 2i3tfi [e 2i( *> - 0) 2 
a ~ 4 1 ^+2 + 2 6 + ^"+2 1 ' 

1,1 2*0 (q + e 2Ut, )e 2iaJt + e 2 ^)e~ 2 ^ |2 

X ~ 4' 2 6 + 2 + a 2 + 2 + /3 2 1 

1, a(a + e 2i4 >)e 2iaJt /3(-/3 + e**) e -WJt 2 

+ 4 1 2 + a 2 + 2 + /3 2 I ' 

1„ 1 (1 + ae**)e* aSt (1 - /3e 2 ^)e" 2i/3j * |2 
2/ = 7(1 - o + + 



4 VI 2 2 + a 2 2 + f3 2 1 

, a(l + ae 2i4l )e 2iajt f3{l - pe**) e -WJt 
+ l 2+^ + 2+^ 1 ( } 

1,1 (1 + ae 2i(t> )e 2iaJt (1 - /3e 2 ^) e - 2i/3j * |2 
~ 4 1 2 + 2W + 2T^ ' ' 
1 J{-(3 + e - 2i ^)e 2iJtf3 a(a + e- 2i *)e- 2iJta 
Z ~ 4^ (3 2 + 2 a 2 + 2 ^ 

, (1 + ae 2i *)e 2iJta (1 - Be 2i *)e- 2iJtfi 1. 
1 a 2 + 2 /3 2 + 2 2 J 

(a + e - 2i *) e - 2iJto + e - 2i *)e 2i3tfi _ 1 _ 2j(/) 

+ ( a 2 + 2 + /3 2 + 2 2 6 } 

/9(1 - p e *<l>) e -xJtfi _ a(l + ae 2i *)e 2iha 
x ( ^T2 ^+2 )} ' 



1 e 2iJt ^{e- 2 ^-/3)P e- 2i3ta a(a + e- 2i4 >) 



X 



/3 2 + 2 a 2 + 2 

,e 2lito (e 2 ^a + 1) e - 2lJ ^{l - e 2i *(3) 1 , 
- a 2 + 2 + /3 2 + 2 2' 



-2a to(a + e - 2 ^ ) i e 2 ^(e- 2 ^-/3) 
+ ( ^T2 2 6 + ^T2 } (A10) 

" 2iJ ^/3(l - e 2 ^/3) e 2iJto a(e 2 ^a + 1) 



x 



2 



/3 2 + 2 a' 
(a + e 2 ^) + l e2 ^ + e- 2 ^(e 2 ^-/3) 



1 a 2 + 2 2 /3 2 + 2 1 

e 2iJto (e 2 ^a + 1) e ~ 2iJtl3 (l - e 2 ^/3) 1 
' ^+2 + ^T2 + 2'^' 
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FIG. 1: (Color online) C(p > A ) vs. and t, for (a) A = 0.9 and (b) A = -0.9. Here the length of 
the chain is chosen to be N = 4. 
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FIG. 2: (Color online) E max as a function of A, for (a) different values of D and (b) corresponding 
0. 
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FIG. 3: (Color online) t op t vs. A, for (a) different values of D and (b) corresponding <fi. 




FIG. 4: Singlet fraction F vs. A at t op t for different values of D. 
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FIG. 5: (Color online) First maximum entanglement vs. D and corresponding <j). Graphs (a) and 
(b) refer to the case A > and (c) and (d) refer to the case A < 0. 
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FIG. 6: (Color online) t op t in terms of D, for different values of (a) A > and (b) A < 0. 
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FIG. 7: (Color online) (a) E max and (b) t op t in terms of inverse temperature for different values of 
D at A = 1. 
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FIG. 8: (Color online) E max and (b) t op t in terms of inverse temperature for different values of D 
at A = -1. 
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FIG. 9: (Color online) v 1 (for an infinite chain) and t op t (for the chain of length N=8) vs. D for 
the case of A = 0.6 and A = —0.6. 
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